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Abstract 

We examine the entanglement creation between two mutuahy independent two-level atoms im- 
mersed in a thermal bath of quantum scalar fields in the presence of a perfectly reflecting plane 
boundary. With the help of the master equation that describes the evolution in time of the atom 
subsystem obtained, in the weak-coupling limit, by tracing over environment (scalar fields) degrees 
of freedom, we find that the presence of the boundary may play a significant role in the entan- 
glement creation in some circumstances and the new parameter, the distance of the atoms from 
the boundary, besides the bath temperature and the separation between the atoms, gives us more 
freedom in manipulating entanglement generation. Remarkably, the final remaining entanglement 
in the equilibrium state is independent of the presence of the boundary. 
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I. INTRODUCTION 



Quantum entanglement has now been recognized as a key resource in quantum information 
science [1], since it plays a primary role in quantum communication [2], quantum teleporta- 
tion [;:>], quantum cryptography [1] and so on. An interesting issue in the discussions for the 
essence of entanglement, which has attracted a lot of attention, is the relationship between 
entanglement and environment. It is known that an environment usually leads to deco- 
herence and noise, which may cause entanglement that might have been created before to 
disappear. However, in certain circumstances, the environment may enhance entanglement 
rather than destroying it [5, 6, 7, 8, 9, 10]. The reason is that an external environment 
can also provide an indirect interaction between otherwise totally uncoupled subsystems 
through correlations that exist. For example, correlations in vacuum fluctuations or fluctu- 
ations at finite temperature can provide such an interaction, when entanglement generation 
is considered in systems in external quantum fields. 

Recently Benatti et al have discussed, in the framework of open systems, entanglement 
generation for two, independent uniformly accelerating two-level atoms interacting with a 
set of scalar fields in vacuum. In the weak coupling limit, the completely positive dynamics 
for the atoms as a subsystem has been derived by tracing over the field degrees of free- 
dom [11], and there it has been shown that the asymptotic equilibrium state of the atoms 
turns out to be entangled even if the initial state is separable. Similar results have been 
obtained by considering two atoms immersed in a thermal bath of scalar particles at a finite 
temperature [12], where, in contrast to Ref. [11], two atoms are assumed to be at a finite 
separation. It is found that for any fixed, finite separation, there always exists a temperature 
below which entanglement generation occurs as soon as time starts to become nonzero and 
for the vanishing separation the entanglement thus generated persists even in the late-time 
asymptotic equilibrium state. Therefore, one can manipulate the entanglement production 
by controlling two controllable parameters: the bath temperature and the separation of the 
atoms. 

In the above studies, the field correlation functions that characterize the fluctuations of 
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fields play a very important role in determining whether entanglement is generated. On the 
other hand, it is well-known that the presence of boundaries in a fiat spacetime modifies 
the fluctuations of quantum fields, and it has been demonstrated that this modification can 
lead to a lot of novel effects, such as the Casimir effect [13], the light-cone fluctuations when 
gravity is quantized [1 1], the Brownian (random) motion of test particles in an electromag- 
netic vacuum [15], and the modification for the radiative properties of uniformly accelerated 
atoms [16]. 

A question then arises naturally as to what happens to the entanglement generation if the 
field correlations are modified by the presence of a reflecting boundary. Now we have one 
more controllable parameter other than the separation and the bath temperature, i.e., the 
distance of the atoms from the boundary and another interesting question is what is the role 
that the new parameter plays in the entanglement generation. These are questions we are 
going to address in the present paper. We shall examine the entanglement generation of two 
non-interacting two-level atoms immersed in a thermal bath of scalar particles subjected to 
a perfectly reflecting plane boundary. With the help of the master equation that describes 
the evolution of the open system (atoms plus external thermal fields) in time, we find that 
the presence of the boundary may play an significant role in controlling the entanglement 
creation in some circumstances and the new parameter, the distance of the atoms from the 
boundary, gives one more freedom in controlling the entanglement generation. It is, however, 
interesting that the probable remaining entanglement for the asymptotic equilibrium state 
at late times is not dependent on the presence of the boundary. 

II. TWO ATOM MASTER EQUATION 

The system we shall examined is composed of two independent two-level atoms in weak 
interaction with a set of massless quantum scalar fields at a finite temperature T. We assume 
that a perfectly reflecting plane boundary for the scalar fields is located at 2; = in space 
and one atom is placed at point xi and the other at X2. Without loss of generality, we take 
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the total Hamiltonian to have the form 



(1) 



Here Hg is the Hamihonian of the two atoms, 



= 2^'^i^ (" = 1,2) 



(2) 



where = (Xj ® ctq, a- = uo ® Cj, CTj, (i = 1, 2, 3) are the Pauh matrices, ctq the 2x2 
unit matrix, n = (ni, 77-2, ^3) a unit vector, uj the energy level spacing, and summation over 
repeated index is implied, if^ is the standard Hamiltonian of massless, free scalar fields, 
details of which is not relevant here and H' is the Hamiltonian that describes the interaction 
between the two atoms with the external scalar fields which is assumed to be weak. The 
general form for H' can be written as 

3 

^' = Y1 [^^/^ ® ^o)$/.(t, xi) + (ao ® X2) ] . (3) 

Now we assume that the scalar fields can be expanded as 

N 

%{^) = E + ix;r<P'^^\x)] , (4) 

a=l 

where 0^^-' (x) are positive and negative energy field operators of the massless scalar field, and 
are complex coefficients that "embed" the field modes into the two-dimensional detector 
Hilbert space and play the role of generalized coupling constants [J I]. It should be pointed 
out that the coupling constant A in (1) is small, and this is consistent with the assumption 
that the interaction of the atom with the scalar fields is weak. 

It is well-known that the evolution of the total system density (i.e., the two atoms plus the 
environment) in time obeys the Liouville equation dtptot{t) = —i[H, ptot{t)] with the initial 
total density having a generic form ptot{0) = p(0)®Pb, where the environment fields are taken 
to be in a thermal state characterized by and the atom in an initial state p(0). Since our 
interest is in the djTiamics for the two atoms only, we must trace over the environment degrees 
of freedom and concentrate on the analysis of the reduced time evolution, p(t) = Tr^[ptot]- 
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Provided that the field correlations decay sufficiently fast at large time separations, or much 
faster than the characteristic evolution time of the subsystem alone, the reduced density of 
the two-atom subsystem can be proven, in the limit of weak-coupling, to obey an equation 
in the Kossakowski-Lindblad form [9, 11, 12, 17] 

dp{t) 



dt 



■[H,f,,p{t)]+C[p{t)] , (5) 



with 



and 



a,f3=l 



1 2 

a,/3=l 

The coefficients of the matrix Cf"'^^ and i^l"^'' are determined by the field correlation func- 
tions in the thermal state p^: 

G:;^(^-o = /3(<^>.(t,x.)<f,(t',x^))^. (8) 

The corresponding Fourier and Hilbert transforms read respectively 

/"OO 

Qt,'\\)= / rfte^^*GSf)(t), (9) 



dtsign{t)e^''G[f\t) = - diu ^^^^ ' (10) 

where P denotes principal value. One can show that the Kossakowski matrix Cl°'^^ can be 
written explicitly as 

cr'= E er'K-)*-!?*,"', (11) 

where 

iplf = UiUj , ip^f = ^{Sij - riiUj ± ieijkrik) ■ (12) 

Similarly, the coefficients of H^j^ can be obtained by replacing QI^'^\^uj) with K}'^i^\e,uj) in 
the above expressions. For the sake of simplicity of our treatment, we now assume that the 
field correlation functions are diagonal such that 

G^it - t') = 5,,G{t - t', X, - x^) . (13) 



This requirement can be fulfilled by demanding that the coupling coefficients satisfy the 
following condition 

N 

E = (14) 

a=l 

or by assuming that the field components ^i{x) are independent. 



III. THE CONDITION FOR ENTANGLEMENT CREATION 



With the basic formalism established, now we shall start to examine whether entangle- 
ment can be generated between two independent atoms in external thermal fields at a finite 
temperature T subjected to a reflecting boundary (i.e. fields are constrained to vanish on 
the boundary), in particular, what is the influence the presence of a boundary that modifies 
the quantum correlations of the fields will have on entanglement generation. 

For the sake of simplicity, let us further assume that two atoms are separated from each 
other by a distance L and are at an equal distance z from the boundary(Fig. (1)), i.e., 
Zi = Z2 = z. Due to the assumption that the fields reflect from the boundary completely, we 
can use the method of images [18] to find the field correlation functions ( Eq. (8) ), 



^ oo 



47r2 



(t — t' — im(3 — ieY 
6i 



{t-f - imf3 - ief - {2zy 



(15) 



Gr(t-t') = GSf(t-t') 



47r2 



^ oo 



_{t - t' - imf3 - iey - 



{t-f - im(3 - ie)2 - (2^)2 - L'^ 
where /? = l/{kT). Plugging Eq. (15) and Eq. (16) into the Eq. (9), we can easily obtain 



(16) 



;(^^)(A) = ^^^^^(A) = ^ _ % ^ sin (22: A) 



2n 1 - e-l^^ 27r 1 - e'l^^ 2zX 



5. 



A 



2n 1 - e-P^ LA 



sin(LA) 5ij A sin( VI^TIi^A) 
2^ 1 - e-/5^ + Az^X ~ 



(17) 
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FIG. 1: The infinite conducting plate is taken to lie along the plane z = 0, the distance between 
each atom and the plate is z, the atom separation is L 



According to Eq. (11), we can write 



a 



(11) 



(12) 



r^{2i) _ 



= AiSij - iBieijkUk + CiTiiTij 
A26ij - iB2eijknk + C2ninj , 



;i8) 



and the corresponding coefficients are 



C2 



LU 1 



4tt 1 



sm(2zu) 



sin(22;u;) 
2zuj 



LU 1 + e 



-fSuj 



47r 1 - e-^" 



sin(Lu;) sin(-\/ + Az'^cu) 



^/iF+Az^u 



Izu) 



Bo 



Air 



sin(Lu;) sin( V-^^ + Az'^uj) 



Luj 



CO 1 + e 
uj 1 + e-^"^ 



1 + 



sm(2zuj) 



Awl 



2zuj 

sin(Lu;) sin(v^L2 + Az^cu) 



Luj 



VL^ + Az^u 



(19) 



7 



Similarly, the /C^ for the Hamiltonian H^s can be obtained easily, but here we do not give 
the formulae in detail. As has already been discussed in detail elsewhere [11, 12], the effective 
Hamiltonian H^s can be expressed as a sum of three pieces. The first two correspond to the 
corrections of the Lamb shift at a finite temperature which should be regularized according 
to the standard procedures in quantum field theory and nevertheless they can be accounted 
for by replacing u in the atom's Hamiltonian Hs with a renormalized energy level spacing 

u = u + t[IC^\-uj) - IC^\uj)] . (20) 

Meanwhile the third is an environment generated direct coupling between the atoms and it 
is temperature independent. So the term associated with H^s in (5) can be ignored, since we 
are interested in the temperature-induced effects. Henceforth, we will only study the effects 
produced by the dissipative part C[p{t)]. 

Using the explicit form of the master equation (5), we can investigate the time evolution 
of the reduced density matrix and figure out whether the state of the two-level atom system 
is an entangled one or not with the help of partial transposition criterion [19]: a two- 
atom state p(t) is entangled at t if and only if the operation of partial transposition of 
p{t) does not preserve its positivity. In general, the two-atom system in the thermal bath 
will be subjected to decoherence and dissipation, which may counteract the entanglement 
production, so that the final equilibrium state is very likely to be separable (however this 
may not always be true as we will demonstrate later). But if we consider the system evolving 
in a finite time, during which the decoherence and dissipation are not dominant, the initial 
separable state may evolve to an entangled one. Here, we adopt a simple strategy for 
ascertaining the entanglement creation at a neighborhood of the initial time t = 0, which 
has been introduced in Ref. [9]. For simplicity, we let the initial pure, separable two-atom 
state be p(0) = |+)(+| ® |— )(— | and consider the quantity 

Qit) = {x\m\x) , (21) 

where the tilde signifies partial transposition and \x) is a properly chosen 4-dimensional 
vector. According to the results of Ref. [9, 12]), entanglement is created at the neighborhood 
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of time t = (i.e., dtQ{0) < 0), if and only if 

{u\C^^^^\u){v\{C^^^Y)\v) < \{u\Re{C'^^^^)\v)\^ 



(22) 



where the subscript T means matrix transposition and the three-dimensional vectors \u) and 
\v) can be chosen in a simple form as Ui = Vi = {1, —i, 0}. Using Eq. (19), we can calculate 
Eq. (22) for the vector n along the third axis directly and deduce that the condition (22) 
becomes 



A2 

A, 



where 



Ai 



-I3lu 



/ sin(Lu;) sin(i^vI?+4? 



/I- 



sin(2ct;z) 



(23) 
(24) 
(25) 



For a given energy gap for the atoms, Bi'^ / Ai'^ takes the values in the interval [0, 1] and is 
only temperature-dependent, while the value of ^42^/^1^ is determined by two parameters, z 
and L and is temperature independent. One can see that when the temperature is zero i.e. 
/3 — *■ CX3 and L is not infinite, the inequality (23) is always satisfied, therefore entanglement is 
generated. At the same time, if the separation is vanishing (L = 0), then / Ai'^ becomes 
unity and the inequality (23) is always obeyed and thus entanglement created too, no matter 
where the atoms are placed, as long as the bath temperature is not infinite. 

Let us now discuss what happens when z ^ ot z L, and z ^ oo or z ^ L, i.e., when 
the atoms are placed very close to and very far from the boundary. Expansion of Eq. (25) 
in power series oi z/L yields 

n 2 
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— uoL cos{ujL) + sin{ujL) 



5cu6L6 



u!Lcos{ujL) — sin(ct;L) 



LuL{-15 + Lu'L') cos(cjL) + 3(5 - 2lu'L') sin(cjL) 



(26) 



The corresponding form for z/ L ^ oo reads 



sin^(co'L) L sin(co'L) 



Z LU 



313 



sm{uL) sm{2uz) — uoL sul{lJilL^/l^!^^z!^^JLP■) 



(27) 



where sm^{ijjL) / uj'^L?' is just the value of A-^ / Ai' without the presence of the boundary [12] 
(i.e., the corresponding value of A2^ / Ai^ in the limit oi z oo). It is interesting to note 
that in the limit of z/L ^ 0, the leading term of A2^ / A^ is independent on z and is only 
a function of ujL (refer to Eq. (26)). This leading term differs from the value of A2^ / Ai' in 
the case without the boundary. As a result, when the atoms are placed very close to the 
boundary, the presence of the boundary will have a significant effect in determining whether 
the inequality (23) is satisfied or whether entanglement is created. In Fig. (2), the leading 
term of A-^ / A^ when z/L ^ is plotted as a function of ujL vs A2^ / Ai^ in the case without 
the boundary. This Figure reveals that when uoL is small, approximately smaller than 3, 
that is when the separation, L, is approximately less than three times the characteristic 
wavelength of the atom's radiation (but L is still large enough to maintain z <^ L), the value 
of 742^/^41^ in the case with the presence of a boundary will be appreciably larger than that 
without as long as L is not vanishingly small. This means that at a certain temperature 
the presence of a boundary would make the atoms be entangled which otherwise still be 
separable. Therefore the presence of the boundary provides us more freedom in controlling 
entanglement creation in this case. However, when ujL is large, i.e., the separation is much 
larger than the characteristic wavelength of the atom's radiation, the value of A2' / Ai' with 
the presence of the boundary generally becomes smaller than that without, since it decreases 
faster (as power of (cuL)"^ as opposed to (cuL)"^) as uoL grows. Therefore, in this case the 
presence of the boundary will make the atoms less likely to be entangled than otherwise. 
Meanwhile when z/L is very large, i.e., when the atoms are very far from the boundary, the 
influence of the presence of the boundary on the entanglement generation is negligible as 
expected and this can be easily seen from Eq. (27) since now the leading term is the same 
as the value of A2^ / Ai^ in the unbounded case. 

To have a better understanding, we also plot, in Fig. (3), A2^ / Ai' as a function of two 
dimensionless variables, z/L and ujL, according to Eq. (25). One can see from this figure 
that, as 2/L varies, appreciable oscillations occur when ujL is of order one and when uoL is 
very small, the value of 742^/^41^ is very close to unity and does not oscillate significantly 
as z/L varies. At same time, A2^ / Ai^ also decays very fast with the increase of ujL and 
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FIG. 2: The dashed hne represents J^j A3[ in the limit of zjL 0, and the solid line denotes 
the function, sin^((jL)/(a;L)^, i.e, A2/A1 without presence of a boundary. Approximately in the 
interval [0,3] of ojL, the value of ^2^/^!^ with the presence of the boundary is always larger than 
that without. 

the oscillations (as z/L varies) is damped dramatically. So we conclude that both when 
uL is very small or very large, the variation of location of the atoms has no significant 
influence on the entanglement generation. Note, however, that this by no means suggests 
that the presence of the boundary does not affect the entanglement generation (refer to the 
discussions in the preceding paragraph). 

The next question we want to ask is what is the maximum difference between the value 
of / Ai'^ with the presence of a boundary and that without. We will try to answer 
the question numerically and approximately. For this purpose, let us plot, in Fig. (4), 
A2^/Ai^ as a function of z/L with a set of fixed values of the dimensionless parameter 
ujL in both the cases with and without a boundary. Our numerical calculations as shown 
illustratively in Fig. (4) indicate that the maximum fluctuation of A2^ / Ai^ in the case with 
the presence of the boundary around that without as a function of z/L is at the neighborhood 
of uL ^ 2.027. Setting cuL ^ 2.027, we find that the function of yl2^/Ai^ vibrates around 
sin^ {cu L) / (cu LY ^ 0.196. Due to the fact that the swing of vibrating function A2'^ / Ai'^ is 
slowly decreasing with the increasing z/L, we can find out the maximum value of A2^/ Ai^, 
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FIG. 3: A'l' I A\' as a function of two dimensionless variables, zjL and ujL. 

by adjusting the parameter zjL^ to be approximately 0.416 achieved at zjL <^ 1. This 
gives the maximum effects of the presence of the boundary on A2^/Ai^ or equivalently on 
entanglement creation. 

To get a more concrete picture, let us take a typical transition frequency of a hydrogen 
atom, u ~ for an example. Then uL ^ 2.027 means L ^ 6.08 x 10"^m which is 

much larger than the usual size of an atom. It is easy to find that in the unbounded space the 
inequality (23) is satisfied or entanglement is created between two atoms if the temperature 
is below 262.663K. However, with the presence of a boundary, we find that the upper bound 
in temperature for entanglement generation can be increased to T < 379.731K. This is a 
hundred Kelvins improvement. Note for L ^ 6.08 x 10~^m, we still have a plenty of room 
to satisfy z/L <ti 1, so the maximum value of 0.416 is used for y42^/Ai^ here. 

Finally, let us briefly discuss what happens if the two-atom system is not aligned strictly 
parallel to the plane boundary. Take the distance from the plane of the atom which is closer 
as z, then the distance of the other atom from the plane will be be larger or smaller than z 
depending on whether the system is inclined away from or towards the boundary. Therefore 
the effect of inclination of the system is that the field correlation function with respect to the 
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FIG. 4: The ratio A'^ j Ax^ is described function of the z/L in the real hne, with the parameter 
are selected a;L = (0.5,1,1.5,2.027,2.8). The dashing line expresses the value of sin^(a;L)/(a;L)^ 
which is the value of ^42^/^1^ without the boundary conditions. As we can easily see that all the 
values of ^42^/^1^ vibrate around the sin^(Lo'L)/(u;L)^. 

atom which is displaced and cross correlation function G\j [t — t') = G\- (t — t') effectively 
get a smaller effective z if the system is inclined towards the plane and a larger effective one 
if otherwise (refer to Eqs. (15,16)). Consequently, taking into account the fact that / Ai^ 
is an oscillating function of z when the system is parallelly placed, one would expect that 
if the two atom system is originally located parallel to the plane at where this function is 
at its peak value the inclination in either direction will make the entanglement creation less 
likely to occur. In contrast, if the system is located at where this function is at its local 
minimum the inclination in either direction will make the entanglement generation more 
likely to happen. However, when the system is placed at any point in the interval where 
function A2^ / Ai^ is monotonically increasing, the atoms will be less likely to entangle if the 
inclination is towards the plane and more likely if otherwise. Similarly, when the system is 
located at any point in the interval where function ^42^/^1^ is monotonically decreasing, the 
entanglement generation will be more likely to come about if the inclination is towards the 
plane and less likely if otherwise. 



13 



IV. THE ENTANGLEMENT OF THE EQUILIBRIUM STATE WITH THE 
BOUNDARY 



In the preceding Section, we find that, in certain circumstances, the presence of a boundary 
plays a significant role in generating entanglement between atoms initially prepared in a 
separable state in a thermal bath of external quantum scalar fields and in fact entanglement 
is created as soon as time starts if the inequality (23) is satisfied. However, the condition (23) 
does not tell us whether the entanglement thus generated can persist in late times, or whether 
the final equilibrium state is still entangled or not. 

At late times, the two-atom subsystem will be in the asymptotic equilibrium state. Though 
the effects of decoherence and dissipation will generically make the state be separable so that 
no entanglement is left in the end, there are also cases in which the entanglement still exists 
at late times. To examine whether the final equilibrium state is entangled or not, let us 
assume, without loss of generality, the reduced density matrix to have the form 



o-Q ® (To + Poi(t) (To ® (Ji + Pio(t)ai ® (To + Pij{t)(Ji ® aj 



(28) 



where the components poi{t), Pioit), Pijit) are real. Substituting Eq. (28) into Eq. (5), we can 
obtain, with setting i/cfr = 0, 

dpQi{t) 



dt 
dpto{t) 



-AA^PQiit) - ABin, - 2B2n,T + 2B2nkP^kit) - 2Cipoi(t) + 2Cin,nkPQk{t) , (29) 
-AAipioit) - ABim - 2B2niT + 2B2nkPki{t) - 2Cip,o(t) + 2CininkPko{t) , (30) 



dt 

dpijit) _ _g^^^^^.^^^ _ 4A2Pji{t) + 4A2T6ij - 4Bi[nipoj{t) + njp,o{t)] 



dt 

-2B2[nipjo{t) + rijpoiit)] + 252K(Pfeo(^) + Pokit))]5ij 
-ACipijit) - AC2Pji{t) + 4:C2[ninkPjk{t) + rijUkPkiit) - riirijr] 
+2Ci[ninkPkj{t) + UjUkpikit)] + 4C2[r - nknipki{t)]5ij . (31) 

Here, r is the trace of the density matrix r = ^^^iPuit). Recall that Ui are the compo- 
nents of the unit vector appearing, for example, in Eq. (2) and Eq. (12). If we symmetrize 
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and anti-symmetrize the density matrix components Poiit), Pijit), we can split the above 
system of differential equations into two independent sets. One can then show that the anti- 
symmetrized components decay exponentially as time grows while the symmetrized ones 
approach a non-zero asymptotic value. Therefore, there exists a final equilibrium state p, 
the explicit form of which, for a non-zero atom separation, can be found by setting the right 
hand side of Eq. (29), Eq. (30), and Eq. (31) to be zero, since any equilibrium state satisfies 
dtp = 0, and the solution is 

i2A, + A2)B^{B,-B2) 

(Ai-A2)5i(2Ai + A2)n, 



POi — PiO 



2Ai' - Ai'A^ - A2B1B2 + Ai{B2' - A 



2\ ' 



. {A^-A2)B^{2B^ + B2)n,n, 

2A^' - Ai^A2- A2B1B2 + Ai{B2^ - A2^)' ^ ^ 

To see if the equilibrium state is entangled or not, we will calculate its concurrence which is 
defined to be C[p\ = max{Ai — A2 — A3 — A4} , where A^ , /i = 1, 2, 3, 4 , are the square roots 
of the non-negative eigenvalues of the matrix pp in decreasing order. Here, the auxiliary 
matrix p = {0-2® 0-2) P^ {0-2® 0-2) ■ Substituting Eq. (32) into Eq. (28), we can easily calculate 
the concurrence for the unit vector n along the third axis via the above the equations. It is 
found that the concurrence is zero, which means that the equilibrium state is separable and 
entanglement generated initially does not persist at late times. Same results can be obtained 
for n along other directions. 

However, it should be pointed out that all the expressions of Eq. (32) become indefinite 
(of the form 0/0) when the separation of the atom approaches zero which results in Ai = 
A2,Bi = i?2,C*i = C*2- Therefore, the case for the vanishing atom separation should be 
dealt with separately. Taking the trace of both sides of Eq. (31) for the vanishing atom 
separation, we find that r = T,f^iPii{t) is actually a constant of motion, which is determined 
by the initial reduced density, while the expression for r in Eq. (32) is no longer valid. In 
fact, the positivity of the initial density matrix requires that —3 < r < 1. Consequently, we 
should take r as a new independent parameter, and components of the density matrix for 
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the equilibrium state, p, in the present case, read 

where R = Bi/Ai. The corresponding concurrence can be calculated directly 



C{p) = max 



5R'-3 

T 



2(3 + L 3 _ /22 
which is non-zero provided r for the initial state p(0) obeys 



, (34) 



This reveals that when the atom separation is zero (L = 0), the entanglement generated 
initially persists at late time despite of the decoherence and dissipation of the external 
environment and the late-time equilibrium state is still entangled, as long as (35) holds. 
This is in sharp contrast with the case of a non-zero separation. However, the presence of 
the boundary has on effect on deciding whether the initially created entanglement can be 
maintained at late times in the equilibrium state, since the concurrence is only dependent 
on r and R, and factors containing the boundary parameter z are all canceled out in the 
expression of R if one recalls Eq. (19), thus the concurrence for the final equilibrium state is 
independent of the presence boundary. 



V. DISCUSSION 



In summary, we have examined the entanglement generation between two mutually in- 
dependent two-level atoms immersed in a thermal bath of scalar particles subjected to a 
perfectly reflecting plane boundary. With the help of the master equation that describes 
the evolution in time of the atom subsystem obtained by tracing over environment (exter- 
nal scalar fields) degrees of freedom, we find that the presence of the boundary may play 
a significant role in controlling the entanglement creation in some circumstances and the 
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new parameter, the distance of the atoms from the boundary, gives one more freedom in 
controlhng the entanglement generation. 

In particular, when two atoms are placed very close to the boundary, i.e., z/L 1 and 
uoL is approximately less than three, that is, when the separation, L, is approximately less 
than three times the characteristic wavelength of the atom's radiation, then for a certain 
temperature the presence of the boundary will make the atoms be entangled which would 
otherwise still be separable. Therefore the presence of the boundary gives us more power 
in creating entanglement. However, when uoL is large, i.e., the separation is much larger 
than the characteristic wavelength of the atom's radiation, the presence of the boundary will 
make the atoms less likely to be entangled than otherwise. Meanwhile when z/L is very 
large, i.e., when the atoms are very far from the boundary, the influence of the presence of 
the boundary on the entanglement generation is negligible as expected. 

At the same time, we find that the variation of location of the atoms has significant 
influence on entanglement generation between two initially independent atom only when uoL 
is of order one, or in different words, both when ujL is very small or very large, the variation 
of location of the atoms has no appreciable effect on the entanglement generation. Note, 
however, that this by no means suggests that the presence of the boundary does not affect 
the entanglement generation. 

Our analysis also reveals that the entanglement generated because of the correlations 
induced by the environment will persist in the late time asymptotic equilibrium state if 
the separation between the atoms is vanishing. However, when the separation is non-zero, 
the entanglement will disappear at late times and the asymptotic equilibrium state becomes 
unentangled again. Finally, the presence of a boundary generally has no effect on maintaining 
the entanglement initially generated in the asymptotic equilibrium state. 
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